arXiv: 1503.07978v2 [cond-mat.stat-mech] 8 Oct 2015 


Tensor product variational formulation applied to pentagonal lattice 


Michal Daniska and Andrej Gendiar^ 

^Institute of Physics, Slovak Academy of Sciences, SK-8f5 11, Bratislava, Slovakia 

The uniform two-dimensional variational tensor product state is applied to the transverse-field 
Ising, XY, and Heisenberg models on a regular hyperbolic lattice surface. The lattice is constructed 
by tessellation of the congruent pentagons with the fixed coordination number being four. As a 
benchmark, the three models are studied on the flat square lattice simultaneously. The mean- 
freld-like universality of the Ising phase transition is observed in full agreement with its classical 
counterpart on the hyperbolic lattice. The tensor product ground state in the thermodynamic 
limit has an exceptional three-parameter solution. The variational ground-state energies of the spin 
models are calculated. 
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I. INTRODUCTION 

The tensor product state (TPS) has been proved to be an appropriate ansatz for obtaining a ground-state of 
strongly correlated quantum systems. Numerous computational approaches and methods have been developed and 
successfully applied to one-dimensional (ID) and two-dimensional (2D) quantum systemJilt^. The main purpose of 
this work is to address the tensor product state to two-dimensional spin systems on a non-trivial lattice configuration, 
which is represented by an infinite hyperbolic surface. The lattice is made from the pentagonal tessellation and forms 
a negatively curved surface of the constant Gaussian curvature. Such a lattice geometry has not yet been considered 
for 2D quantum systems in the thermodynamic limit. Lacking information on constructing an appropriate numerical 
algorithm for the pentagonal lattice is given. 

The adherence to the mean-field-like universality class was concluded in the 2D classical Ising systems on various 
types of hyperbolic lattices. The mean-field-like behavior of the phase transition exponents is caused by the non- 
Euclidean geometry of the underlying lattice^^^^, and is not related to the analysis by the mean-field approximation. 
The classical spin systems have been analyzed by the Corner Transfer Matrix Renormalization Group (CTMRG) 
method^, which is a variant of the DMRG for ID quantum system^. Now, we extend our earlier studies, which dealt 
with the classical spin systems, to quantum spin systems on the curved pentagonal lattice surface. 

We, therefore, propose an approximative sche me of the round-a-face TPS, which is related to the standard vertex 
type TPS with two-state auxiliary variableJ^^. Our main intention is to approximate the ground state of spin-| 
models by a minimal number of variational parameters. Improvements of the numerical accuracy go beyond the scope 
of this work. 

The paper is organized in the following way. Section II specifies the model Hamiltonians on the Euclidean and 
the non-Euclidean lattices. The variational approach is discussed with respect to the TPS approximation. A short 
description of the numerical algorithm and discussion of tensor symmetries are included within two subsections. The 
numerical results are analyzed in Sec. Ill, and we comment the results in Sec. IV. 


II. THE MODEL 

We calculate the ground-state energy of the quantum Ising, XY, and Heisenberg models on two different lattices: 
the Euclidean square lattice and the hyperbolic pentagonal lattice. Figure [^depicts the two lattice types constructed 
by tessellation of the congruent polygons: either by the squares {p = 4) shown on the left or by the pentagons {p = 5) 
on the right. The alternating gray-white color scheme of the polygons is chosen for improving the visibility only. The 
spin variables are located on the vertices of the regular polygons with the constant coordination number, which is 
equal to four in both cases. We study properties of the spin system in the thermodynamic limit, i.e., the number of 
the lattice vertices is infinite. The Hamiltonian H of the three models defined on the both lattices can be expressed 
in the following compact form 


n = ( 1 ) 

It), 

(v) 

where Gjf represents the local Hamiltonian of the p-sided polygon, the lattice is constructed from, and k marks the 
position of the polygon on the lattice. The summation runs over all the positions of the polygons (fc)^. The polygon 
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Figure 1: The Euclidean square lattice geometry {p = 4) and the pentagonal hyperbolic lattice geometry (p = 5). Notice that 
the hyperbolic lattice is made of the regular pentagons of the equal size and shape. Nevertheless, the mapping of the lattice 
onto the Poincare disk depicts them deformed with exponentially decreasing sizes toward the circle boundary. 



Figure 2: Graphical representation of the local Hamiltonian with its particular shape; the square (p = 4) on the left and 
the pentagon (p = 5) on the right. 


on the position is described by the ordered set of spin indices fci, ^ 2 , kp, see Fig. where ki stands for the 
unique number which is assigned to the corresponding vertex within the labeling scheme of the lattice vertices. The 
local Hamiltonian has the expression 
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where , S^,, are the Pauli operators, and the spin indices obey the cyclic condition kp+i = ki. The x- 
component of the external magnetic field is described by the variable h, and the spin couplings J^y and Jz specify 
the three models. In particular, J^y = —Jz = 1, and h = 0 describe the Heisenberg model, Jxy = 1 and Jz = h = 0 
specify the XY model, whereas J^y = 0, = 1 with an arbitrary h lead to the transverse field Ising modePliSI. Since 

the ferromagnetic ordering of the models leads to a simpler TPS formulation, we have opted for the positive coupling, 
Jxy = 1, and we consider dimensionless units throughout the entire work. The constant prefactors ^ and ^ reflect 
the sharing of the spin couplings and the magnetic field, respectively, if the Hamiltonian is formed by the polygonal 
tessellation in Eq. 0. 

Our objective is to obtain the ground-state of the system 


|$p) = lim 

^ N^oo 
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in the thermodynamic limit by a variational minimization of the ground-state energy normalized per bond 


c(p) 
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where N stands for the total number of the lattice spins, aj, j = 1, marks one of the two base states I or 

t of the lattice spin and N}, denotes the total number of the bonds (the nearest-neighbor pairs). In order to 
simplify the numerical calculation, we approximate |$p) by a TPS iT'p), which is given by the product of the identical 
tensors Wp of the same polygonal structure as either of the two local Hamiltonians has (cf. Fig.0. The p-rank 
tensors depend on p spin-4 variables labeled by indices ki, ...,kp with two base states ctj,. = j, or f. Tne p individual 
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spin variables are grouped into a single one with 2^ base configurations denoted as {cr^} to simplify the notations 
if necessary. It means that the tensor element Wp{{ak}) = bbp(o’fci'^fe 2 " ' )■ For instance, there are 32 base spin 

configurations for the pentagons, which can be represented in the arrow notation as l iiiiil, {'tiiti}) •■•j 

{ttttt}. Thus, the approximative ground state in the form of the polygonal has the following form in the 

thermodynamic limit 


<7i(T2‘‘‘<^N (^)p 


(5) 


where the sum runs over the 2^ base spin states. Since the TPS |'I'p) has the product structure of the identical 
tensors Wp, the variational problem in Eq. is in the thermodynamic limit equivalent to the minimization of the 
local energy of an arbitrary bond in the lattice center (to avoid boundary effects) 
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where £ is the index of a polygon containing the selected central bond and the factor 2/p reflects that each polygon 
contains p bonds shared with neighboring polygons. Moreover, the product structure of Idtp) enables us to express 
the denominator 


{%\%) = E n %*(K})<J{<}.{^jW^p({a4) = V{Wp{{a})) 
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and the numerator 
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as sole functions of the tensor elements Wp({cr}), where we removed the subscript k due to the uniform TPS. Here, 
stands for the corresponding matrix element of the local Hamiltonian G^^\ is the Kronecker 

symbol, and (fc)^ \ {£} denotes the set of all polygon indices except for the index £. 

Consequently, the minimization over the set of variational parameters in Eq. is replaced by a much 

simpler problem 


= 


2 U{Wp{{a})) 
WAW}) p V{Wp{{<7})) 


(9) 


where we minimize over 2^ tensor elements Wp{{a}) only. This set can be further reduced if additional constraints 
are taken into account. The translational and rotational symmetries shrink them to three variational parameters, 
provided that the ground state jd/p) is evaluated at zero magnetic field, i.e., no symmetry-breaking mechanism is 
present in the system. The calculation of the numerator N{Wp{{a})) and the denominator V{Wp{{a})) in Eq. Q is 
carried out separately by means of a numerical algorithm described below. 

Thus defined TPS approach offers a new perspective of estimating the 2D quantum systems on hyperbolic surfaces. 
Such a significant approximation is sufficient only for the case of p = 5 since we have observed a weak entanglement 
and non-critical behavior of the classical Ising model on variety of hyperbolic lattices in our earlier studies®^]. An 
exponential decay of the density matrix spectra and the correlation function result in the phase transition, where the 
correlation length ^ 1 is finite, reaching its maximal value at the phase transitiorP^. Eor this reason, we apply the 

current scheme of TPS as it can be a sufficient approximation for the ground-state properties of the models on the 
pe ntag onal lattice. (However, the square lattice requires to consider a larger dimension in the tensors W 4 , cf. Refs.^ 
ancPT) 


A. The TPFV algorithm 


We refer to the algorithm as the Tensor Product Variational Eormulation (TPVF) and it consists of two parts. 
The first one evaluates the ratio in Eq. 0 by applying the CTMRG methocP separately to the numerator and 
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Figure 3: The double-layer tensor structure of on the left and ^sda^crj.}) on the right. 


the denominator for a given set of the tensor elements Wp{{a}). The second part contains a multi-dimensional 
minimizer, the Nelder-Mead simplex algorithrrP^lli^, which uses the first part to search for the optimized set of the 
tensor elements Wp({o’}), which minimize the ratio in Eq. (|^. The minimizer starts from an initial simplex in the 
space of free variational parameters, one vertex of which is specified by the initial tensor elements Wp{{(j}). The 
simplex undergoes an iterative sequence of size changes and moves towards lower energies and stops if the energy in 
Eq. (|^ converged. 

The CTMRG was originally developed to study 2D classical spin systems. Primarily, it evaluates the partition 
function (preferably, in the thermodynamic limit) as the configuration sum of the tensor product of the Boltzmann 
weights Wb ■ The central idea of this study is to apply CTMRG to quantum systems by replacing the concept of the 
Boltzmann weight Wb from the classical statistical mechanics by the tensors Wp. In order to do this, let us introduce 
a double-layer tensor Zp with the tensor elements 

Zp{Wk^k}) ^ w/;({4})^{.'},{.,}%({^fc}). (10) 

Notice that there are 2"^^ double-layer base spin configurations {c^.cr^}. Figure [^graphically depicts the double-layer 
tensors Zp at the position k, where each shaded polygonal area represents Wp. Thus, in the language of the classical 
statistical mechanics, the general expression for the denominator I?(Wp({(T})) = (ffplTp) in Eq. Q corresponds to a 
tensor product object, which is equivalent to the partition function of a (non-physical) classical Hamiltonian given by 
the product of the tensors Zp. 

This generalization of the Boltzmann weight enters the CTMRG algorithm for both the lattice geometries and the 
consequent numerical calculation yields the denominator 'D{Wp{{a})) for the given set of the tensor elements Wp{{<7}) 
according to Eq. 0. A similar approach can be also used to determine the numerator M{Wp{{a})), as it differs from 
'D{Wp{{a})) only by the additional double-layer structure at the central position £ containing the local Hamiltonian 
. The detailed survey of the CTMRG algorithm is explained in Refs.*^*^ for both the square and the pentagonal 
lattices. We point out only the most relevant part of the CTMRG algorithm in the following. 

The concept of CTMRG is built on corner transfer tensors and transfer tensors Tj^\ where j = 1, 2, ...,p. Each 
of the (corner) transfer tensors is composed of the polygon representing tensors Zp via recurrence relations indexe d by 
the iteration step n (as specified later). Further details of the construction of the transfer tensors are given in Refs.SEl. 
As a result, the (corner) transfer tensors represent a specific lattice sector formed by the corresponding polygons. 
The square and the pentagonal lattices are constructed from a central polygon surrounded by the alternating sectors 
represented by the tensors and 7)^"^ The central polygon is represented by the tensor Zp or WpG^/'^ 


Wp in 


the structure of the denominator 'D{Wp{{a})) and the numerator NjWp{{a\)), respectively. Figure H illustrates the 
situation for the case of 'D{Wp{{a})). Consequently, the relations Q and ® for the denominator VlWpda})) and 
the numerator N{Wp{{a})), if formulated in the CTMRG language of the (rorner) transfer tensors, take the form 


’^{WpiW})) = lim Tr 



( 11 ) 


MiWM"})) = Tr U"GPw, 

The tensors and 7)^"^ expand their sizes iteratively, following the recurrence relation^^ 
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Figure 4: (Color online) The graphical illustration of the denominator in Eq. (Ill for p = A (left) and p = 5 (right). The central 
polygon tensor Zp is surrounded by p identical (corner) transfer tensors and 7)^"^ bordered by the red thick curves within 
the iteration step in accord with Eqs. (131 and (141. The light blue and the white polygons, respectively, correspond to 
and with the gradually decreasing iteration step (n) towards the lattice boundary, where n = 1. Note that the tensors 

and are not localized on single polygons only (as briefly depicted), but all of the tensors for n = 1, 2,..., (n — 1) including 

Zp are recursively inherited in accord with Eqs. (131 and (14|. 


where the transfer tensors satisfy the cyclic condition = 7^*'"^- The infinite TPS geometry is built up gradually 
by increasing the iteration step n = 2, ...oo, which induces an exponential increase of the degrees of freedom of 
the (corner) transfer tensors. The calculations are kept numerically feasible by means of the renormalization group 
proce dure , which integrates out the least probable configurations in the tensors determined by the reduced density 
matrijl^. 

The recurrence relations are initialized in the first iteration step, n = 1 by construction of the tensors and 

-y-f"—1) tensor Zp for all j = 1, 2, ...,p. If the spin variables ak^ are explicitly included in the tensors, we have 


^5({CTfcCrfc}), 

= ^4({0-feCrj,}) , 

<y'i <7. 

K4 k^ 

«5 ^5 

= 2^4({crfcCrfc}) . 


(15) 


B. The tensor symmetries 

As mentioned earlier, let the polygon base spin configuration {cr^} = {'^ki^k^ '' ’^kp) given in the spin-arrow 
notation. Let function Oj return all such configurations for which the integer number j counts the number of the 
spins aligned upward as listed in Tabs. |T|and[TT| Each line in the Tables contains such spin configurations, which are 
identical with respect to the rotational symmetry operations of the p-sided polygon. We make a difference between 
the spin configurations Oj and 0jt when j = 2 (or j = 3 for the case of p = 5) as they are not rotationally equivalent. 
The prime symbol is used for such configuration (and its rotations), in which either the two (j = 2) or three (j = 3) 
spins, aligned upward in each polygon, cannot be grouped together (or equivalently, if the alternating spin alignment 
is maximal). 

As a consequence, we count six or eight distinguishable configurations do • • • for p = 4 and p = 5, respectively. 
These configurations determine the number of the variational parameters Wp{9j) used in the calculation of the 
transverse field Ising model. If, however, the spontaneous symmetry-breaking does not affect the solution, the total 
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number of the variational parameters decreases down to four parameters for the Hamiltonian models (for both lattices 
p). Therefore, the spin-inversion symmetry reduces the variational parameters due to the absence of the preferred 
direction of the spontaneous magnetization. In particular, the spin configuration probability for the pairs 9q and Op, 
9i and Op-i, etc. becomes identical, which can be formally generalized into the equations 

— «min{j,p-j}’ 

^ ^ (16) 

WpiOp) = WpiO^p.py) = 1. 

There are the three equations in the upper expression defining the new variational parameters Oq, Qi, 02 , and one 
equation in the lower expression for 02 /, which has already been eliminated from the set of the free parameters by 
putting 02 / = 1 being the normalization condition in Wp. Hence, the three free variational parameters suffice to 
approximate the ground-state wave function of the models with no spontaneous symmetry-breaking phases. For 
the same reason, if we consider the system without the spin-inversion symmetry, there are either five or seven free 
variational parameters for p = 4 or p = 5, respectively, excluding the normalization parameter Wp(02') = 1- 

The number of the operations performed by the implementation of the Nelder-Mead algorithm we used scales 
linearly with the number of the free variational parameterJi^. Of course, increasing the number of the parameters 
prolongs the computational time, and may encounter numerical instability caused by trapping the system in a local 
minimum of the energy, rather than approaching to the correct global minimum, which corresponds to On 

the other hand, a faster Nelder-Mead optimization ena bles to improve the accuracy by increasing the number of the 
CTMRG states kept (the effective block spin state^^^ and make TPVF more efficient. 


III. NUMERICAL RESULTS 

As we have mentioned earlier, we consider two types of the lattices. The Euclidean one is defined on the regular 
square lattice characterized by p = 4, which serves as a reference for the hyperbolic pentagonal lattice with p = 5 (cf. 
Fig. 0 , which is of our main interest. First, we start analyzing the ground-state properties of the transverse field Ising 
model (TFIM) in details, where the number of the variational parameters Wp{0^) is either five (for p = 4) or seven 


Table I: The full spin configurations for p = 4 sorted by 6j. The last column lists the reduced parameter set with I = 
min{j,p-j}. 
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Table II: The spin configurations for p = 5 as in Tab|T] 
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Figure 5: (Color online) The magnetic field dependence of the six variational parameters IV 4 {9j) on the reference square lattice 
(p = 4). The singular behavior corresponds to the phase transition critical field ~ 3.158. 


(for p = 5), respectively, excluding the normalization parameter Wp(02') = 1- We analyze the TFIM and its phase 
transition by the following three ways: (1) by evaluating the optimized parameters Wp{6j), (2) by the expectation 
value of the magnetization constructed from the ground state jlfp), and (3) by the ground-state energy per bond 

and its second derivative. Second, we reduce the number of the free variational parameters down to three (if all 
symmetries are considered) and calculate the ground-state energies of the XY and Heisenberg models. 

In our numerical calculations we keep at most 20 effective block spin stated, which is sufficient to reach more than 
six-digit convergence in A/’(Wp({cr})) and 'D{Wp({a})). Notice that the numerical accuracy is mainly given by the 
uniform TPS approximation. Further improvement of the numerical accuracy requires to consider a non-uniform TPS 
(i.e. violation of the translational lattice symmetry), and a gradual expansion of the tensor order by implementation 
of auxiliary variable^^^ is necessary. This is, however, beyond the scope of our interest. 

The Nelder-Mead algorithm has been tested for various initial conditions, most of them leading to the identical 
results. Typically, we started with a simplex, where we set all the coordinates of one of its vertices to unity (the 
vertex coordinates represent the free variational parameters Wp{0j)). The initial simplex size was set to O.P^ESl. 

Figures and illustrate behavior of the optimized variational parameters Wp{6j) in the TFIM with respect to 
the externm magnetic field h for the particular lattice geometry p. It is evident that above a certain magnetic field 
(depicted by the vertical dotted lines), specific pairs of the variational parameters collapse onto identical values, 
which yields four single curves in this region for both the lattices. These pairs are exactly those coupled by the 
system of equations in Eq. (16), i.e. those representing spin configurations which are equivalent with respect to spin 
inversion. The magnetic field, at which the collapse causes a singular behavior of Wp(0 ), corresponds to the quantum 


phase transition of TFIM at the magnetic field which we analyze below. In the ordered phase at h < hc^', the 
distinct optimized values of the coupled parameters Wp{9j), as specified in Tabs. [T| and [Hj reflect the existence of the 

spontaneous symmetry-breaking in the TFIM for both the lattice types. In the disordered phase at h> h^\ the four- 
parameter description coincides with the free variational parameters Og, a^, 02 , and the normalization parameter 02 /. 
This confirms the relevance of the additional symmetries in such systems, where the spontaneous symmetry-breaking 
mechanism is not present, such as in the XY and Heisenberg systems at the zero magnetic field. 


(p) 



















Figure 6 : (Color online) The magnetic field dependence of the eight variational parameters W 5 {6j) on the hyperbolic pentagonal 
lattice [p = 5) with the phase transition field ~ 3.264. 


Having evaluated the optimized free variational parameters, we can easily reconstruct the approximative eigenstate 
I flip) as the TPS and apply it for the evaluation of the spontaneous magnetization 


( 5 “) 




(17) 


where a = x or z. The expectation value of the spin operator S^. is evaluated in the central part of the lattice in 
order to suppress all boundary effects. Here, (Sp) denotes the order parameter of TFIM and specifies the quantum 
phase transition at the phase transition field. The resulting dependence of magnetization (Sp) and (Sp) with respect 
to the magnetic field h is shown in Fig. The quantum phase transition is reflected by the singular behavior of all 
magnetization profiles. 

Analyzing the TFIM, we calculated = 3.158 for the Euclidean square geometry and = 3.264 for the 
hyperbolic pentagonal lattice. The most relevant value of the critical magnetic field for the TFIM on the 2D Euclidean 
lattice (as obtained by a recent Tensor RG algorithm) yields = 3.043!^^^. Since the TPS is built up by the tensors 
W 4 of the too low dimension, the long-range order correlations are excludecP, and the mean-field-like dependence near 
the phase transition is necessarily observed. The expectation value of the magnetization {Sp) obeys the scaling relation 


{s;{h)) ^ - hy\ ( 18 ) 

The mean-field exponent /3p = ^ is observed regardless of the lattice geometry. As a reference, the numerical TRG 
analysis gives correct = 0.3295 on the square lattice (free of any mean-field approximations), which is also in 

agreement with Monte Carlo simulation^l^. The inset of Fig. displays the magnetization squared, which confirms 
the mean-field-like behavior for both p by its linear dependence on h if approaching the phase transition. The blue 
dashed lines serve as guides for the eye in order to enhance the linearity. 

A more detailed analysis of the TPVF approximation near the phase transition can be visualized by evaluating the 
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Figure 7: (Color online) The transversal {Sp) and the parallel (Sp) magnetizations with respect to the magnetic field for both 
lattice geometries. The inset shows the mean-field dependence of the magnetization when approaching the transition field; the 
linearity is depicted by the blue dashed lines. 


effective (field dependent) exponent which converges to fip when approaching the phase transition field h[ 


(p) 


I3p= lini /3^^^(/i) = lini 


dHS;{h)) 

a In -h) 


(19) 


as plotted in Fig. The effective exponent obviously converges to the mean-field exponent Pp = \ for both lattice 

types if the phase transition field is approached from the ordered phase, i.e., if ln(/ii^^ — /i) —>■ —oo. The inset shows 
the same dependence on larger scales. The critical exponent on the square lattice (the black curve for p = 4) starts 
deviating at around h > 2.0 from the expected exponent (estimated by the blue dashed curve), which is known to 
converge to = 0.329^^^. 

The ground-state energies normalized per bond given by Eq. © are depicted in Fig. There is a tiny 
difference in both of the energies. Therefore, we plotted the energy difference in the inset (which has no 

physical meaning) in order to point out the difference. The two vertical lines in the inset correspond to the phase 
transitions (for reference only). The magnetic susceptibility 


Xp = 


02 4 ^) 

dh^ 


( 20 ) 


is an independent physical quantity, which characterizes the phase transition. Figure shows the functional depen¬ 
dence of the susceptibility on the magnetic field. The calculation of Xp requires a precise data for the ground-state 
energy because of performing the second derivative numerically. The shape of the non-diverging discontinuity in Xp 

(p) 

at the phase transition magnetic field he is the typical consequence of the mean-field universality. 

The XY and the Heisenberg models are critical (at p = 4) at zero magnetic field. Therefore, the full tensor 
symmetries are considered, cf. Tabs. |l] and |TT[ resulting in the three fr ee v ariational parameters Oq, a^, and 02 , which 


is sufficient to obtain the ground-state energy per bond Eg \ Table 


III 


variational tensor parameters and the energies. The comparison of Eg 
results by SandvilP^*^^ serves as an estimate of the exact ground-state energy £g^^ for the XY and the Heisenberg 
models on the hyperbolic pentagonal lattice, which has not been determined yet. 


summarizes thus calculated optimized free 

on the square lattice with the Monte Carlo 
(5) 
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Figure 8: (Color online) The detailed dependence of the inverse effective magnetic exponent on the magnetic field in the 
logarithmic form. The inset shows the effective exponent at wider magnetic field scale. The blue dashed line estimates behavior 
of the correct effective exponent for the Euclidean lattice. 


Table III: List of the optimized free variational parameters aj which minimize the ground-state energy per bond on the 
square and pentagonal lattices for t he XY and Heisenberg models. The ground-state energies on the square lattice are 
obtained by Monte Carlo simulation^^^^. 


lattice 

model 

p = 4 

p = 5 

XY 

Heisenberg 

XY 

Heisenberg 

“0 

0.69333 

0.50746 

0.66001 

0.48553 

0.1 

0.88010 

0.74826 

0.83844 

0.71804 

“2 

0.94236 

0.83422 

0.90733 

0.80263 

TpiP) 

-1.0846 

-1.3089 

-1.0815 

-1.2913 

ttiMC 

i^Q 

-1.09765 

-1.33887 

— 

— 


IV. CONCLUSIONS 

We have applied the TPVF algorithm to the calculation of the quantum systems (transverse field Ising, XY, and 
Heisenberg models) on the pentagonal hyperbolic lattice with the constant negative Gaussian curvature and the 
coordination number hxed to four. As the reference lattice, the Euclidean square lattice was taken. We have used 
the RG-based algorithm, which approximates the ground state as the product of identical tensors representing the 
congruent polygons; either squares {p = 4) or pentagons {p = 5) for the Euclidean or hyperbolic lattices, respectively. 
We have analyzed the symmetries in the tensors. The uniform TPS reduces the infinite number of the variational 
parameters in the thermodynamic limit down to six (for p = 4) or eight (for p = 5) for the Ising model in the 
presence of the transverse magnetic held. If no spontaneous symmetry-breaking mechanism is present, the number 
of the free variational parameters approximating the tensor product ground state further shrinks down being three 
only. We have analyzed the phase transition of the Ising model by (I) the numerical calculation of the optimized 
variational parameters, (2) the expectation value of the spin polarization, (3) the phase transition magnetic exponent, 
(4) ground-state energy and the magnetic susceptibility. The mean-held-like universality class was observed in the 
vicinity of the phase transition magnetic helds . The ground-state energy per bond has been evaluated for 
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Figure 9: (Color online) The ground-state energy of the TFIM with respect to the magnetic field h. The small difference 
between the two energies is plotted in the inset. 



h 

Figure 10: (Color online) The magnetic susceptibility Xp of fh® TFIM as a function of h for both lattices. 


the XY model and the Heisenberg model for p = 4 and p = 5. 

Our earlier studies of the classical spin lattice models on various types of the hyperbolic surface^l^^ exhibited a 
couple of interesting features, which are important to mention here. The classical Ising model on the pentagonal 
lattice results in a higher phase transition temperature if compared to the Euclidean square lattic^. Analogously, the 
TFIM model in this work also exhibits the identical feature, i.e., Classical spin systems on any hyperbolic 

lattice belong to the mean-field-like universality class in accord with the current study. The mean-field-like behavior 
of the exponents in the classical systems originates from the hyperbolic lattice geometry, not from the numerical 
CTMRG method; the CTMRG accurately reproduces the critical exponents on the 2D Euclidean lattices xhe 
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low dimension of the tensors in the TPS approximation of the round-a-face type suppresses the quantum long-range 
correlations on the square lattice {p = 4) near the criticality, where we obtained h^c'^ deviating by 3.7% if compared to 
the recent numerical analysiJi^. The discrepancy is the consequence of the low dimension of the tensor W 4 appearing 
in the ground-state approximation, which limits the numerical accuracy if approaching the criticality, and it results 
in the mean-field exponent /3 = |. Therefore, it is essential to distinguish between the mean-field approximation of 
the method used and the mean-field-like behavior caused by the hyperbolic lattice geometry. That means that we 
have dealt with both of the ’’mean-field” aspects in the current study. The physical reasoning of the mean-held-like 
behavior caused by the hyperbolic lattice comes from the exceedance of the critical dimensionality (being dc = 3 or 
dc = 4, respectively, for the classical or quantum systems with the nearest-neighbor spin couplings). Notice that the 
Hausdorff dimension of the pentagonal lattice is infinite. 

Strictly speaking, classical spin systems exhibit the phase transitions exclusively in the center of the infinite hy¬ 
perbolic lattices because of strong boundary effects. The ratio between the number of the spins on the boundary 
with respect to the number of the spin on the remaining lattice area inside is « 2.73 for the pentagonal lattice. For 
this reason, the hyperbolic lattice becomes significantly sensitive to the boundary effect^^. On the other hand, the 
boundary effects on the Euclidean lattices play no role in the thermodynamic limit when the lattice size is expanded 
to infinity. If calculating the correlation length ^ on the pentagonal lattice (the distance between two spins measured 
along a geodesics) for the classical Ising model, no divergence at the phase transition was present, i.e., 0 < ^ < 1^. 
Moreover, the correlation function was found to decay exponentially even at the phase transition, which is the direct 
consequence of the hniteness of the correlation length at the phase transition (and is related to an analogous expo¬ 
nential decay of the density matrix spectra)P. In that sense, the phase transition on pentagonal hyperbolic surface is 
non-critical, and we use the term critical for the spin models on the square (Euclidean) lattice only since it is related 
to the divergence of the correlation length by definition. The identical non-critical behavior for p = 5 has also been 
confirmed in this work (for instance, by evaluating the non-divergent magnetic susceptibility at It is highly 

non-trivial to evaluate the correlation length unambiguously as both of the ’’mean-field” aspects are closely related 
and almost impossible to be clearly separated or distinguished. 

An exponentiall y fas t decay of the reduced density matrix spectra has been observed for the classical systems on 
hyperbolic surface^^^. Similar behavior for the quantum systems had to be taken into a possible scenario in the 
current work. For all the above-mentioned reasons, we conjecture the TPFV analysis of the models on the hyperbolic 
pentagonal lattice is more accurate than on the Euclidean ones. It is caused by the mean-held-like behavior (the mean- 
held-like universality) as the consequence of the pentagonal hyperbolic lattice with the inhnite Hausdorff dimension, 
which has been studied by the numerical improved mean-held approximation. Our preliminary results of the identical 
quantum Hamiltonian on different hyperbolic lattices (to be published elsewhere) support our claims and are in 
complete agreement with the conjectures we have made for the classical spin systems. 
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